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This  report  has  been  generated  under  the  advanced  gun  propellant 
evaluation  analysis  portion  of  Project  62602F  2560  and  is  presented  as  a  § 

state  of  the  art  advancement  that  can  be  applied  to  gun  propellant  research 
efforts.  The  study  was  conducted  during  the  period  February  1970  to 
January  1971. 

This  technical  report  has  been  reviewed  and  is  approved. 


CHARLES  PETRI  DBS 

Chief,  Advanced  Development  Division 
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ABSTRACT 


computer  program  is  presented  for  computing  the  chemical  equilibrium 
reaction  products  ^sociated  with  gun  propellant  combustirr.  This  program 
wil!  provide  a  good  first  approximation  of  flame  temperature,  specific 
heat  ratio,  impetus,  and  mean  mo.’ jcular  weight  of  the  combustion  products, 
he  program  ^  developed  at  the  Air  Force  Armament  Laboratory  was  based  on 
a  code  ^ed  by  the  National  Aeronautics  and  Space  Administration  Lewis 

changes  to  the  program  are  discussed,  some 
results  from  advanced  solid  propellant  formulations  calculations  are  pre¬ 
sented,  and  a  form  of  users  manual  for  thermochemical  programs  and  for  the 
program  to  generate  specific  heat  polynomials  from  JANAF  tabular  data  is 
presented. 
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SECTION  I 


INTRODUCTION 


The  field  of  gun  propellant  development,  which  had  been  relatively 
static  for  many  years,  has  recently  begun  to  attract  increased  attention. 
The  availability  of  an  accurate  theoretical  thermochemical  propellant 
performance  computer  program  is  of  central  importance  to  informed 
research  and  development  efforts  in  this  field.  As  a  formulation  screen¬ 
ing  tool,  a  few  dollars  of  machine  time  can  save  hundreds  of  dollars  of 
laboratory  small-batch-mix  effort  by  eliminating  unpromising  candidate 
formulations  and  providing  a  good  first  approximation  to  the  primary 
important  propellant  parameters.  These  would  include  flame  temperature, 
specific  heat  ratio,  impetus,  and  mean  molecular  weight  of  the  combustion 
products .  To  most  expediently  obtain  a  computer  program  for  computing 
the  chemical  equilibrium  reaction  products,  the  Air  Force  Armament 
Laboratory  modified  an  existing  free-energy  minimization  rocket  propellant 
performance  program.  The  program  selected  was  the  one  used  by  the 
National  Aeronautics  and  Space  Administration  Lewis  Research  Center.  This 
code^^i  has  a  background  of  almost  ten  years  of  successful  application. 

In  addition,  it  was  readily  available  and  the  JANAF  specific  heat  data 
could  be  easily  applied  to  the  thermodynamic  data  tape.  The  program,  as 
modified  for  gun  propellant  applications,  was  checked  against  established 
experimental  data  and  other  analytic  chemical  equilibrium  computing 
techniques,  with  excellent  agreement.  Additionally,  an  algorithm  was 
added  which  computes  combustion  gas  viscosity  and  thermal  conductivity 
by  a  Leonard-Jones  potential  technique. 

The  following  sections  of  this  report  discuss  the  minor  changes 
made  to  the  Lewis  approach,  give  some  results  from  advanced  solid  propel¬ 
lant  formulation  calculations,  and  provide  a  form  of  users  manual  for  the 
therraochemical  programs  and  for  the  program  to  generate  specific  heat 
polynomials  form  JANAF  tabular  data^^^ .  Appendix  I  presents  the  Lewis 
program  listing  and  Appendix  II  presents  the  polynomial  fit  program  used 
at  the  Armament  Laboratory. 
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SECTION  II 

FREE  ENTHALPY  EQUATIONS 


This  report  will  not  discuss  the  mathematical  techniques  of  chemical 
equilibrium  computations.  This  task  has  effectively  been  accomplished  in 
a  recently  issued  monograph^^J  which  discusses  in  detail  not  only  the 
Lewis  technique  but  also  several  other  commonly  used  approaches.  This 
report  will  demonstrate  only  that,  at  chemical  equilibrium,  the  Gibbs 
free  energy,  G,  (or  fre*  enthalpy)  will  be  a  minimum. 

From  the  First  Law  of  Thermodynamics,  in  an  irreversible  process: 

dE  +  pdv  -  Tds<0  (1) 

or  at  equilibrium 

dE  +  pdv  -  Tds=0  (2) 

For  a  constant  pressure  and  temperature,  equation  (2)  becomes 

d(E  +  pv  -  Ts)  =0  (3) 


=  H  -  T,S, 
T  T  11 


Equation  (4)  establishes  the  relation  that  the  Gibbs  free  energy  will 
be  a  minimum  at  equilibrium,  and  when  coupled  with  the  atomic  species  mass 
conservation  equations,  it  provides  a  sufficient  set  for  solution  of  the 
combu  tion  species  concentration  and  flame  temperature  (if  the  combustion 
pressure  is  specified) .  Then,  to  set  up  a  minimization  algorthm  only 
requires  a  value  for  the  heat  of  formation  of  the  reactant  compounds  and 
tabular  or  polynomial  data  for  the  specified  heat,  C  enthalpy,  Hj,  and 
entropy,  S-p,  of  the  possible  combustion  products.  These  thermodynamic 
data  can  be  obtained  from  reference  2  for  most  combustion  products  of 
interest:  as 


»T  -  Hjss  =  Cj,  dT 

(5) 

^  dT 

(6) 
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l- 


4' 
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Selecting  a  polynomial  form  for  C^,  the  terras  and  can  be 
conveniently  fitted  to  similar  polynomials  by  performing  the  indicated 
integrations 

-^  =  aj  +  a2T  +  33?^  +  ^47^  + 

^  =  a,  +  1/2  32?  +  1/3  33?^  +  1/4  347^  +  1/5  agT^  +  ag/T  (8) 

g 

^  =  ajLnT  +  3.^1  +  1/2  ajT^  +  1/3  34!^  +  1/4  agl'^  +  ay  (9) 

For  a  more  precise  fit,  the  Lewis  program  utilizes  two  sets  of 
coefficients  for  each  species,  one  in  the  temperature  range  0“  to  1000°K 
and  the  other  from  1000“  to  6000“K.  Limits  of  species  existence  can  be 
placed  on  the  propellant  data  to  restrict  the  range  of  polynomial  appli¬ 
cability. 

The  primary  change  made  to  the  Lewis  algorithm  was  to  additionally 
develop  an  approximate  isochoric  temperature  from  the  isobaric  temper¬ 
ature  provided  by  the  program.  This  requires  the  assumption  that  approx¬ 
imately  the  scime  species  would  form  either  in  an  isochoric  or  isobaric 
process,  providing  a  similar  total  heat  release.  If  these  assumptions 
are  warranted,  then: 


“  S’^Fp 

(10) 

Qr  =  CyTpv 

(11) 

Tpv  =  t'p/CyTpp  =  YTpp 

(12) 

Once  the  isochoric  flame  temperature  and  mean  molecular  weight  of 
the  combustion  products  are  known,  the  impetus  is  simply 


where  R  =  2780  ft-lbs/lb-mole-‘’K 


'^FV 

M 


(13) 


The  program  also  computes  the  constitutive  products  mole  fractions 
as  a  function  of  the  pressure  ratio  (or,  for  a  gun  system,  of  the 
velocity)  if  quasi-isentropic  conditions  are  assumed. 

Equation  (13)  indicates  the  importance  of  developing  advanced  pro¬ 
pellants  with  lower  molecular  weight  combustion  products,  if  the  impetus 


of  a  propellant  is  to  be  increased  without  a  corresponding  flame  temper¬ 
ature  increase.  High  propellant  flame  temperatui  s  are,  of  course,  the 
prime  limiting  factor  in  developing  an  acceptable  barrel  life  for  auto¬ 
matic  aircraft  cannon  having  a  high  performance  and  a  high  rate  of  fire. 
Current  single-base  solid  gun  propellants  display  molecular  weights  in 
the  order  of  25.  Table  I  illustrates  the  distinctive  flame  temperature 
advantage  of  lower  molecular  weights  in  obtaining  cooler  combustion  con¬ 
ditions  without  sacrificing  energy  output. 

TABLE  I.  PLANE  TEMPERATURES  AS  A  FUNCTION  OF 
MOLECULAR  WEIGHT  AND  IMPETUS 


Flame  Temperature,  °K,  at  Impetus,  Ft-lb/lb 


Molecular  Weight, 

Atomic  Mass  Units  325,000  350,000  375,000  400,000 


The  very  favorable  performance  data  indicated  by  the  lower  right- 
hand  entries  in  Table  I  are  responsible  for  the  continued  advanced  solid 
gun  propellant  development  efforts.  The  test  cases  described  in  the 
following  sections  have  resulted  from  calculations  of  two  of  the  formu¬ 
lations  now  under  contract  investigation  and  from  a  comparison  of  results 
of  two  theoretical  computations  >  ith  experimental  measurements  of 
inservice  propellants.  All  computations  are  performed  at  a  standard 
reference  pressure  of  5,000  psia. 

The  first  requirement  for  free  energy  computations  is  a  reliable 
list  of  heats  of  formations  for  the  constituents  to  be  considered.  Table 
II  provides  this  data  for  the  ingredients  in  the  formulations  discussed 
in  this  report.  Additional  tables  are  available  in  references  4  and  5. 
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SECTION  III 


COMPin’ATION  OF  VISCOSITY  AND  CONDUCTIVITY 

For  gun  propellant  applications,  it  is  often  quite  important  to  determine 
a  theoretical  value  for  the  combustion  gas  conductivity  and  viscosity. 
These  are  particularly  necessary  if  the  convective  gun  barrel  heating 
characteristics  of  the  propellants  are  to  be  evaluated.  To  accomplish 
this  task,  the  Leonard-Jones  potential  techniquef^)  was  applied. 

The  viscosity  for  a  single  constitutive  species  is  given  by: 

=  2.6693  X  10‘*  /aHl  (14) 


where  T  is  the  temperature,  M  is  the  molecular  weight,  and  a  and  are 
the  Leonard-Jones  force  constants.  The  viscosity  of  the  mixture  is 
defined  by 

*  M  K  m  4 


with  n  being  the  number  of  species  considered  and  x^  being  the  mole 
fraction  of  the  species.  The  value,  is  defined  as 

The  thermal  conductivity  of  a  species  is  given  as: 


K;  =  1.9891  X  10' 


The  mixture  value  is  given  as 


X  . 


J 


(18) 


This  approach  is  strictly  valid  for  only  dilute  mixtures  of  non-polar 
molecules.  Five  species  are  considered,  (one  of  which,  H2O,  is  polar), 
and  all  have  a  dense  gas  system.  Thus,  the  results  contain  only  order  of 
magnitude  validity,  but  are  none  the  less  of  value,  considering  the  state 
of  the  art  in  mathematical  modeling  of  the  transient  heat  convection 
process  in  gun  bores.  Table  III  gives  the  species  considered  which  are 
typically  98  to  99  percent  of  the  combustion  gas  cunstitutients . 
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TABLE  III.  LEONARD- JONES  Pi 

JTENTIAL  PARAJiETERS  (FORCE  CONSTANTS) 

Species 

Molecular 

Viscosity, 

Conductivity, 

Weight, 

Atomic  Mass 
Units 

Grams/cm-sec 

Cal/cm-sec-K 

H. 

2.016 

2.915 

38.0 

28.02 

3.749 

79.8 

CO 

28.01 

3,.  706 

88. 

CO2 

44.01 

3.897 

213. 

H2O 

18.0 

2.824 

230.9 

Table  IV  consists  of  the  fi-JOl  relations  needed  for  the  solution  of 
Equations  (14)  and  (17).  The  firsc  column  contains  values  for  the 
relatively  non-polar  species  such  as  H2,  N2,  CO,  and  CO2,  while  the 
second  column  has  a  corrected  value  for  the  highly  polar  H2O. 
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TABLE  IV. 

—  -  -  —  —  - _ i 

REDUCED  CROSS  SECTION 

KT/  e 

NON -POLAR 

POLAR 

KT/e 

— 

NON -POLAR 

POLAR 

1.0 

1.5938 

2.6199 

7 

0.8725 

0 . 8946 

1.2 

1.4568 

2.4257 

8 

0.8536 

0.8422 

1.4 

1.3557 

2.2713 

9 

0.8378 

0.8043 

1.6 

1.2800 

2.1413 

10 

0.8242 

0.7760 

1.8 

1.2216 

2.0263 

11 

0.8123 

0.7544 

2.0 

1.1751 

1.9217 

12 

0.8017 

0.7376 

2.2 

1.1377 

1.8258 

13 

0.7922 

0.7243 

2.4 

1.1066 

1 . 7373 

14 

0.7836 

0.7134 

2.6 

1.0803 

1.6559 

15 

0.7756 

0.7045 

2.8 

1.0579 

1.5812 

16 

0.7883 

0.6971 

3.0 

1.0385 

1.5126 

32 

0.6939 

0.6462 

3.2 

1.0214 

1.4499 

64 

0.6262 

0.6033 

3.4 

1.0063 

1.3926 

128 

0.5634 

0.5528 

3.6 

0.9928 

1.3401 

256 

0.5056 

0.5006 

3.8 

0.9807 

1.2922 

512 

0.4528 

0.4505 

4.0 

0.9696 

1.2485 

5.0 

0.9265 

1.0797 

6  0 

0 . 8960 

0.9694 
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SECTION  IV 


INPUT  FOR  GUN  PROPELLANT  PROGPJVM 


The  inpjJt  data  for  the  gun  propellant  program  (Appendix  I)  can  be  divided 
into  four  general  groups  hr/ing  code  names  as  follows: 

1)  Thermodynamic  data  for  the  reaction  products  (THERM(3  data) . 

2)  Data  pertaining  to  the  reactants  or  propellants  (REACTANTS  data) . 

3)  Special  options  related  to  chemical  specios  present  in  the  combusted 
gas  (0MIT/INSERr  data)  . 

4)  Namelist  data  including  the  type  of  problem,  pressure  ratio  and  area 
ratio  schedules,  various  options,  etc.  (NAMELISTS  data). 

The  required  order  of  the  data  cards  is: 

1)  One  card  with  the  code  word  THERM0  punched  in  columns  1  to  6. 

2)  THERM0  data. 

3)  REACTANTS  card.  Number  of  reactants  right  adjusted  in  columns  1  to 
3  and  number  of  cases  right  adjusted  in  columns  4  to  6. 

4)  0MIT  and/or  INSERT  data. 

5)  One  card  containing  up  to  80  columns  of  alphanumeric  identification 
data. 

6)  NAMELISTS  data. 

For  any  particular  problem  there  may  be  multiple  REACTANTS  cards.  Each 
set  of  REACTANTS  cards  may  be  followed  by  multiple  NAMELISTS  input  cards. 
Each  type  of  input  data  is  discussed  in  detail  in  the  following  sections. 

1 .  'niERHjJ  DATA 

The  thermodynamic  data  for  the  reaction  products  may  be  read  either  from 
cards  or  from  magnetic  tape.  IVhen  tape  input  is  i;sed,  both  the  THERM0 
code  card  and  all  THERM0  data  cards  must  be  omitted.  It  is  anticipated 
that  most  users  will  prefer  the  magnetic  tape  since  this  reduces  the 
number  of  cards  which  must  be  handled.  The  card  format  for  the  THERM0 
data  is  described  in  Section  V.  Use  of  the  tape  input  requires  no  action 
on  the  part  of  the  user,  once  the  required  data  is  established  on  the 
tape. 

2.  REACTANTS  DATA 

REACTANTS  data  is  required  for  all  problems.  Following  the  code  card 
there  should  follow  one  Ccird  for  each  reactant  species  being  considered, 
with  a  maximum  of  IS  allowed.  Each  reactant  card,  after  the  code  card, 
contains  the  following  information: 


a.  The  cliemical  formula  for  the  species. 

b.  Either  the  number  of  moles  of  the  reactant  or  the  relative  weight. 

The  relative  weight  is  the  weight  of  each  fuel  or  oxidizer  expressed 
as  a  fraction  of  the  total  fuel  or  oxidizer. 

c.  The  enthalpy  of  the  species  expressed  in  calories  per  mole.  (This  is 
not  required  for  an  assigned  temperature  problem.  See  Section  IV. 4 
for  description  of  assigned  temperature  problem.) 

d.  The  state  of  the  reactant  (gas,  liquid,  or  solid). 

e.  The  temperature  associated  with  the  enthalpy. 

f.  Designation  of  each  reactant  species  as  either  a  fuel  or  an  oxidizer. 
(The  program  then  combines  all  fuels  into  an  effective  fuel  and  all 
oxidizers  into  an  effective  oxidizer.) 

g.  The  reactant  species  density.  (This  information  is  optional  and  may 
be  omitted.  If  densities  are  input  for  all  reactants,  the  program 
calculates  a  density  for  the  overall  system.) 

h.  Multiple  cases  with  identical  reactants  are  identified  by  the  reactants 
code  card  and  the  preceding  data  need  not  be  repeated. 

The  reactant  information  is  arranged  on  a  card.  Standard  chemical 
alphabetic  symbols  are  required  for  the  chemical  elements  (i.e.,  H  for 
hydrogen,  HE  for  helium,  LI  for  lithium,  BE  for  beryllium,  etc.).  Each 
reactant  may  be  composed  of  no  more  than  five  distinct  chemical  elements. 
Each  chemical  element  is  allowed  two  columns  (colmms  1  and  2,  10  and  11, 

19  and  20,  28  and  29,  and  37  and  38)  for  its  symbol.  For  those  elements 
whose  symbol  consists  of  a  single  letter,  the  symbol  must  be  placed  in  the 
left-most  (left  adjusted)  of  the  two  columns.  The  formula  number  (the 
number  of  atoms  of  each  element  in  a  reactant)  is  to  appear  in  the  columns 
immediately  following  the  chemical  symbol  (i.e.,  columns  3  to  9,  12  to  18, 
21  to  27,  30  to  36,  and  39  to  45).  The  formula  niomber  may  appear  anywhere 
within  the  seven  column  field,  and  a  decimal  point  must  be  used  for  each 
nmnber.  Note  that  the  exponent  numeric  form,  or  E  format,  may  not  be  used 
here.  Also  note  that  for  a  given  reactant  an  individual  element  may 
appear  only  once  (empirical  formula).  Thus,  a  compound  such  as  furfuryl 
alcohol,  whose  formula  might  be  written  as  C4H30’CH20H  to  represent  its 
structure,  must  be  treated  2is  though  it  were  C5H602  or  02H6C5,  etc.  The 
order  of  the  elements  is  immaterial  for  this  program. 

The  relative  weight  of  each  reactant  should  be  placed  in  columns  46  to  52. 
The  number  may  appear  anywhere  in  these  columns,  a  decimal  point  must  be 
used,  and  the  exponent  numeric  form  may  not  be  used.  If  the  number  of 
moles  0.'  each  species  is  input  in  these  columns  the  same  restrictions 
apply,  but  an  additional  bit  of  information,  an  M  in  column  53,  must 
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appear  on  the  card.  Column  53  is  left  blank  for  relative  weights.  Note 
that  the  relative  weight  and  molar  designations  may  not  be  mixed  for  any 
single  problem. 

The  reactant  enthalpy  is  placed  in  columns  54  to  62,  and  the  units  must 
be  calories  per  mole.  The  number  may  appear  anywhere  in  these  coluirins, 
a  deci.iial  point  must  be  used,  and  the  exponent  numeric  form  may  not  be 
used.  A  sign  is  required  only  if  the  enthalpy  is  negative.  Obviously, 
this  enthalpy  must  be  consistent  with  the  enthalpy  base  used  for  the 
THERli0  data.  For  the  JANAF  data,  the  enthalpy  base  is  as  follows: 

For  eacii  element,  the  phase,  condensed  or  liquid,  which  is  most  stable 
at  one  atmosphere  pressure  is  designated  as  a  reference,  or  ba'^e,  state. 
The  heat  of  formation  of  this  state  at  298.15®K  is  arbitrarily  specified 
as  zero.  Then  all  thermodynamic  properties  for  other  phases  of  the 
element,  as  well  as  for  any  compound  containing  the  element,  are  refer¬ 
enced  to  the  base  state. 

The  initial  state  for  each  reactant  is  indicated  by  an  S  for  solid,  an  L 
for  liquid,  or  a  G  for  gas  in  column  63.  'fhe  initial  reactant  temper¬ 
ature,  in  degrees  Kelvin,  goes  in  columns  64  to  71.  Both  the  initial 
state  and  temperature  are  used  only  to  label  the  output  printout  and  if 
they  are  omitted  or  input  incorrectly,  the  calculations  are  not  affected. 

For  all  problem  types  (see  Section  IV. 4)  except  a  detonation  problem, 
columns  73  to  70  are  reserved  for  the  reactant  density  in  grams  per  cubic 
centimeter.  The  inputting  of  this  data  is  optional,  but  if  it  is  given, 
a  decimal  point  must  be  used.  For  a  detonation  problem,  columns  73  to  80 
are  used  for  the  reactant  specific  heat  at  constant  pressures.  This  data 
is  required  for  each  reactant  for  a  detonation  problem 

3.  0MIT/ INSERT  DATA 

The  purpose  of  the  0MIT/INSERT  option  is  as  follows.  The  program 
considers  as  possible  products  all  species  in  the  THERM0  data  v/hich  are 
consistent  with  the  cliemical  system  which  was  input  on  the  REACTANTS 
cards.  Occasionally,  it  may  be  desired  to  eliminate  one  or  more  species 
from  consideration  as  possible  products.  This  may  be  accomplished  by 
placing  the  species  chemical  formula  on  an  0MIT  card.  Both  gaseous  and 
condensed  species  may  be  omitted. 

Only  the  ricses  of  cc'ndensed  species  may  appear  on  an  INSERT  card.  The 
INSERl’  option  has  been  included  for  two  reasons.  The  less  important  of 
the  two  is  that  if  it  is  known  beforehand  that  one  or  more  particular 
condensed  species  will  be  among  the  final  equilibrium  composition  for  the 
first  assigned  point  (i.e.,  the  combustion  diamber  for  a  rocket  problem), 
then  a  small  amount  of  computer  time  can  be  saved  by  using  an  IN.  •  RT  card. 
The  more  important  reason  for  the  INSERT  option  is  that,  in  rare 
instances,  it  is  impossible  to  obtain  convergence  without  the  use  of  an 
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INSERT  card.  This  occurs  when,  by  considering  gases  only,  the  temperature 
becomes  extremely  low.  When  this  happens,  the  program  outputs  the  follow¬ 
ing  message:  "LOW  TEMPERATURE  IMPLIES  CONDENSED  SPECIES  SHOULD  HAVE  BEEN 
INCLUDED  ON  AN  INSERT  CARD"  and  no  solution  will  be  obtai.'ed.  The  user 
should  resubmit  the  problem  with  one  or  more  appropriate  condensed  species 
listed  on  an  INSERT  card. 

The  use  of  0MIT  and  INSERT  cards  is  completely  optional;  however  when 
used,  each  card  must  contain  either  the  code  word  0MIT  in  columns  1  to  4 
or  INSERT  in  columns  1  to  6  and  the  chemical  formulas  of  from  one  to  four 
product  species.  The  names  of  the  species  must  be  left-adjusted  and  begin 
in  columns  16,  31,  46,  or  61.  In  addition,  the  species  name  must  be 
punched  on  the  card  exactly  as  it  appears  in  the  THERM0  data.  This  means 
that  the  order  gf  the  chemical  elements  within  the  species  name,  the 
formula  numbers,  and  the  state  of  the  species  must  all  be  identical.  For 
example,  gaseous  ammonia  would  commonly  appear  as  NH3.  Because  of  the 
alphabetizing  scheme  used  for  thermodynamic  data,  this  species  is  written 
as  H3N1(G) . 

4.  NAMELISTS  DATA 

The  word  NAMELIST  is  a  F0RTRAN  IV  statement  and  provides  a  convenient  way 
0^  inputting  data.  One  of  the  major  advantages  of  the  NAMELIST  input 
method  is  that  the  data  can  be  placed  in  any  order  and  need  not  be  placed 
within  pre-specified  card  columns  as  is  the  case  for  the  REACTANTS  data. 

A  secondary  advai.'tage  is  that  the  variable  name  appears  with  its 
numerical  value  on  a  card,  thus  making  it  easy  to  tell  the  purpose  of  each 
piece  of  data  without  having  to  refer  to  a  program  manual  or  to  the 
program  listing.  Even  though  inputting  data  by  means  of  NAMELIST  is 
convenient,  the  following  rules  must  be  observed. 

All  data  input  by  means  of  a  NAMELIST  statement  must  be  associated  with 
a  NAMELIST  name,  and  for  this  program  there  are  only  two:  INPT2  aii. 
RKTINP.  The  first  column  of  all  NAMELIST  data  cards  may  not  be  used  and 
must  be  blank.  The  second  column  of  the  first  card  only  miast  contain  a 
$  character.  The  applicable  NAMELIST  name  must  be  placed  on  the  first 
card  only.  The  name  must  begin  in  column  3,  and  must  be  followed  by  at 
least  one  blank  column. 

Each  NAMELIST  name  is  assigned  a  group  of  F0RTRAN  variables.  (The 
variables  belonging  to  INPT2  and  RKTINP  are  listed  in  Tables  V  and  VI, 
respectively.)  These  variables  a- e  assigned  to  a  NAMELIST  name  within 
the  program,  and  numeric  data  cai,  r-ly  be  input  in  conjunction  with  a 
NAMELIST  name. 

Although,  in  general,  there  are  seven  forms  that  input  data  may  take, 
only  three  are  of  interest  for  this  report.  These  are  integer  constants. 
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real  number  constants,  and  logical  constants.  These  are  defined  as 
fol lows : 

a.  An  integer  consists  of  one  to  eleven  digits  written  without  a  decimal 
point . 

b.  A  real  number  constant  may  be  written  in  two  ways.  The  first  form  is 
from  one  to  nine  digits,  including  a  decimal  point,  followed  by  an 
exponent.  The  exponent  is  written  as  the  letter  E  followed  by  a  sign 
and  a  one  or  two  digit  integer.  For  example  12345.  could  be  written 
as  1.234SE04  or  .12345+5,  etc.  The  second  form  is  one  to  nine  digits, 
including  a  decimal  point,  but  not  followed  by  an  exponent. 

c.  A  logical  constant  may  be  either  true  or  false.  There  are  two  forms 
in  which  the  constant  may  be  written.  These  are  either  .TRUE,  or  T 
.FALSE,  or  F.  The  periods  before  and  after  the  long  forms  are 
mandatory. 

Logical  constants  may  be  associated  only  with  logical  variables.  Integer 
constants  may  Le  associated  with  real  variables  and  vice  versa,  and  the 
poper  conversion  is  made  automatically  by  the  computer.  Note  that  when 
inputting  data  without  a  NAMELIST  statement  this  conversion  may  cause 
difficulties  and  should  only  be  used  with  great  care. 

On  a  NAMELIST  card,  -each  piece  of  data  must  be  written  in  conjunction 
with  its  variable  name  (e.g.,  P  =  100 ,’  or  INDEX  =  3,  or  PSIA  =  .TRUE., 
etc.).  The  equal  sign  is  mandatory,  and  the  data  items  must  be  separated 
by  commas.  If  more  than  one  card  is  required,  the  last  item  on  each  card, 
except  for  the  last  card,  must  be  a  constant  followed  by  a  comma.  A 
variable  name  may  not  be  the  last  item  on  the  card.  The  end  of  a  group 
of  data  is  signaled  by  a  $  character.  This  may  be  on  the  same  card  as 
that  containing  the  NAMELIST  name  if  only  one  card  is  used  or  may  be  on 
any  succeeding  card,  but  it  may  not  be  the  first  character  c.n  a  succeeding 
card.  Note  that  it  is  not  necessary  to  completely  fill  a  card  with  data 
before  beginning  a  subsequent  card.  In  fact,  the  user  may  want  to  put 
each  piece  of  data  on  separate  cards.  This  facilitates  subsequent  data 
changes.  One  fi-.  1  note,  it  is  not  necessary  to  input  data  for  every 
variable  contained  in  a  NAMELIST;  however,  each  variable,  .as  been  assigned 
a  value  within  the  program,  and  this  value  is  used  unless  superseded  by 
input  data.  The  assigned  values  for  all  NAMELIST  variables  in  this 
program  are  included  in  Tables  V  and  VI. 

To  illurtrate  the  preceding  instructions,  the  following  paragraphs  discuss 
the  specific  NAMELIST  data  required  by  this  program,  what  each  data  bit 
does,  and  then  some  specific  examples. 
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The  INPT2  NAMELIST  data  must  be  used  for  all  problems.  The  variables  in 
the  INPT2  NAMELIST  are  listed  and  defined  in  Table  V  along  with  the  type 
and  assigned  value  of  each  variable.  The  type  of  problem  (one  of  either 
TP,  HP ,  RKT ,  or  DETN) ,  at  least  one  pressure  (P) ,  and  the  relative  amounts 
of  fuel  and  oxidizer  (i.e.,  one  of  either  EQRAT,  OF,  FPCT,  or  FA)  are 
required  for  each  problem.  The  other  variables  are  either  completely 
optional  or  may  be  required  only  if  a  particular  problem  type  is 
requested. 

The  variable  JANF  is  a  flag  to  indicate  which  thermodynamic  data 
polynomial  should  be  used.  JANF=0  designates  a  polynomial  of  the  form 

(y=Ai)-A2T+A2T"^+\4T^+A5T^.  This  is  the  form  used  by  the  Lewis  P.esearch 
Center  and  is  also  the  form  the  Armament  Laboratory  uses  after  converting 
the  JANAF  data. 

The  Lewis  program  is  quite  flexible  in  the  type  of  problems  it  will  con¬ 
sider  or  the  thermodynamic  variables  which  may  be  fixed.  For  gun  pro¬ 
pellant  problems,  it  is  most  advantageous  to  assign  the  enthalpy  and 
pressure  which,  in  system  nomenclature,  constitutes  an  RKT  problem. 

a.  RKT  Problem 

If  RKT  =  .TRUE,  or  T,  then  a  rocket  problem  is  solved.  A  rocket  problem 
is  one  for  which  the  combustion  takes  place  at  an  assigned  enthalpy 
(REACTANTS  cards)  and  assigned  pressure  (single  value  in  the  P  variable) , 
followed  by  an  isentropic  expansion.  The  combustion  temperature,  the 
thermodynamic  properties  and  rocket  performance  are  calculated  for  the 
chamber,  the  throat,  and  assigned  exit  points.  The  exit  points  desired 
are  specified  in  the  PKTINP  namelist.  If  a  temperatqre  is  assigned 
(T  variable) ,  thep  combustion  will  be  assumed  to  occur  at  that  temperature 
and  not  at  the  assigned  enthalpy.  The  RKT  problem  type  combines  the  H,S 
and  T,S  problems  of  the  earlier  program  version.  Note  that,  when  this 
option  is  selected,  the  second  set  of  NAMELIST  data,  RKTINP,  is  then 
required.  The  combustion  gas  conductivity  and  viscosity  are  computed  for 
this  option  only. 

b.  DETN  Problem 

If  DETN  =  .TRUE,  or  T,  then  the  program  solves  a  Oiapman-Jouget  detonation 
problem.  The  thermodynamic  properties  and  composition  downstream  of  the 
detonation  wave  are  calculated  for  an  assigned  pressure  and  temperature 
upstream  of  the  wave.  This  detonation  calculation  is  primarily  oriented 
for  gas  phase  detonations.  It  is  of  limited  Utility  for  computation  of 
detonations  in  solids,  as  the  specific  heat  polynomial  of  the  reactants 
must  be  on  the  data  tape,  which  will  not,  in  general,  be  the  case  for  a 
t>pical  solid  explosive  constituent  ingredient. 
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c.  P,  PSIA,  and  MMHG  Variables 

The  value  P  may  be  either  the  combustion  pressure  of  a  rocket  engine 
(RKT  problem)  or  a  pressure  at  which  composition  and  thermodynamic 
properties  (TP  or  HP  problem)  are  desired.  Up  to  26  values  may  be  input, 
but  only  one  is  permitted  for  an  RKT  problem.  The  program  assumes  that 
the  units  of  the  pressures  are  in  atmospheres  unless  either  PSIA  or  MMIG 
are  set  true.  If  PSIA  =  .TRUE,  or  T  then  the  pressure  units  are  assumed 
to  be  in  psia  units.  If  MNHG  =  .TRUE,  or  T  then  the  pressure  units  are 
assumed  to  be  in  millimeters  of  mercury.  PSIA  and  MMIG  may  not  both  be 
set  true  in  the  same  problem. 

d.  T  Variable 

Values  for  the  T  variable  should  be  input  only  if  TP  is  set  true  or  if  an 
assigned  chamber  temperature  is  desired  for  RKT  problem.  The  temperature 
units  must  be  degrees  Kelvin. 

e.  EQRAT,  0F,  FPCT,  and  FA  Variables 

The  relative  amounts  of  total  fuel  and  total  oxidizer  may  be  expressed  in 
one  (and  only  one  for  each  problem)  of  four  ways.  The  mixture  ratio  C0F) 
and  percent  fuel  (FPCT)  are  self  evident.  The  fuel-to-air  weight  ratio 
(FA)  and  equivalence  ratio  (EQRAT)  are  new  to  this  version  of  the  program 
and  are  not  of  use  for  gun  propellant  calculations. 

f.  I DEBUG  Variable 

If  IDEBUG  =  .TRUE,  or  T  then  a  printout  of  the  intermediate  calculation 
details  is  obtained.  This  is  intended  primarily  as  a  debugging  aid  in  the 
event  there  are  problems.  The  output  is  extensive  and  is  explained  on 
page  34  of  Reference  1. 

g.  I0NS  Variable 

If  I0NS  =  .TRUE,  or  T  then  the  program  will  consider  ionized  species. 

The  earlier  Lewis  program  version  could  not  handle  ionic  species  although 
provision  was  made  in  one  of  the  subroutines  for  future  consideration  of 
ions.  This  capability  is  included  in  the  newer  program  and  thermodynamic 
data  for  ionized  species  are  available  as  well  as  a  sample  problem  whi-'-h 
uses  I0NS  =  T. 

h.  EQL,  FR0Z,  PCP,  SUBAR,  and  SUPAR  Variables 

The  five  RKTINP  variables  (EQL,  FR0Z,  PCP,  SUBAR,  and  SUPAR)  in  the  RKTINP 
namelist  are  listed  in  Table  VI.  Of  these,  only  PCP  is  required.  The 
RKTINP  namelist  is  required  only  if  RKT  is  set  true  in  the  INPT2  namelist. 
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(1)  PCP  Variable 

The  value  PCP  is  the  ratio  of  the  chamber  pressure  to  the  exit  pressure. 
This  is  one  of  the  independent  variables  (along  with  the  entropy)  for  the 
nozzle  portion  of  a  rocket  problem,  and  as  such,  at  least  one  value  must 
be  input,  ^s  many  as  22  values  may  be  input;  however,  unlike  the  older 
program,  values  for  the  chamber  and  throat  should  not  be  input.  Pressure 
ratios  for  the  subsonic  portion  of  the  nozzle  may  be  used.  Values  must 
be  ordered  so  that  the  magnitudes  increase  monotonically . 

(2)  SUBAR  and  SUPAR  Variables 

Gun  propellant  performance  at  specific  area  ratios  corresponding  to  a 
given  gas  velocity,  either  subsonic  or  supersonic,  may  be  requested. 
Subsonic  ratios  are  read  as  SUBAR  and  supersonic  as  SUPAR.  IVhen  assigned 
area  ratios  are  requested,  the  range  of  PCP  values  should  be  large  enough 
to  include  the  assigned  area  ratios.  If  the  range  is  not  large  enough, 
extrapolation  may  be  performed.  SUBAR  values  should  be  ordered  such  that 
the  values  decrease  and  the  SUPAR  data  should  increase.  Thirteen  of  each 
variable  are  permitted.  Note  that  the  use  of  this  program  feature  is 
completely  optional. 

(3)  EQL  and  FR0Z  Variables 

The  program  will  calculate  performance  for  both  equilibrium  and  frozen 
expansions  unless  instructed  otherwise.  If  EQL  =  .FALSE,  or  F  then  the 
equilibrium  expansion  is  omitted  and  only  frozen  performance  is  calcu¬ 
lated.  If  FR0Z  =  .FALSE,  or  F  then  the  frozen  expansion  is  omitted  and 
only  equilibrium  performance  is  calculated. 


TABLn  V.  VARIABLES  IN  INPT2  NAMELIST 


Variable  Name 

Variable  Type 

Value  (Unless  Set 
Differently) 

1 

Definitions 

1 

JANF 

Integer 

1. 

Thermodynamic  data 
polynomial  indicator 

P 

Real 

0. 

Assigned  pressures 

T 

Real 

0, 

Assigned  temperatures 

EQRAT 

Real 

0. 

Equivalence  ratio 

0F  1 

1 

Real 

0. 

Oxident-to-fuel  weight 
ratio 

FPCT 

Real 

0. 

Percent  fuel  by  weight 

FA 

Real 

0. 

Fuel-to-air  weight  ratio 

TP 

Logical 

False 

Assigned  temperature  and 
pressure  problem 

HP 

Logical 

False 

Assigned  enthalpy  and 
pressure  problem 

RKT 

Logical 

False 

Rocket  problem 

DETN 

Logi  cal 

False 

Detonation  problem 

PSIA 

Logical 

False 

Assigned  pressure-t-o  be 
in  psia  units 

M^tlG 

Logical 

False 

Assigned  pressures  in 
millimeter  of  mercury 
units 

I0NS 

Logical 

1 

False 

Consider  ionic  species 

I  DEBUG 

Logical 

False 

Print  intermediate 
output 

1 

ABLE  VI.  VARIABLES  IN  RKTINP  N. 
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Table  VIl  contains  the  input  data  for  several  propellant  performance 
cases,  and  the  output  is  listed  on  Tables  VIII  through  XI.  Thesf*  cases 
are  respectively  M-10  (Case  100),  M-9  (Case  101),  an  HMX- rubber  propellant 
system  (Case  102)  ,  and  an  RDX-triamino-guanidine  nitrate  formulation  (Case 
106) .  A  comparison  of  these  data  with  experimental  or  independent  theore¬ 
tical  solutions  gives  the  following  results. 

Case  100 


Tv 

Fp 

Y 

Mw 

Lewis 

3312 

366,548 

1.22 

25.1 

Experimental 

3034 

346,180 

1.23 

•  24.2 

Case 

101 

Lewis 

3853 

394,180 

1.17 

27.2 

Experimental 

3840 

396,840 

1.20 

26.5 

Case 

102 

Lewis 

2566 

377,584 

1.28 

18.9 

Other 

Theoretical 

2545 

373,000 

(-) 

18.9 

Case 

106 

Lewis 

2637 

390,870 

1.29 

18.8 

Other 

Theoretical 

2647 

391,000 

(-) 

18.8 

The  input  data  (Table  VII)  are  discussed  in  the  program  writeup  except 
that  no  0MIT  or  INSERT  cards  are  used. 
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(Tlie  reverse  of  this  page  is  blank) 
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TABLE  XIII.  THERM0  DATA  OUTPUT  FORMAT  FOR  TITANIUM  DIOXIDE 


TITANIUH  OIOXIOE  (TI02>  IDEAL  CAS 
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F— .WW3S6€>Aa-^. 
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900 r-  iMfr—i 
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.95760211E-ai 
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1.130  71.S0(  71.509 


.120469676>02  >.S0104624C-06  .0 

.22037632E>01  .107734792*00 
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TABLE  XIV.  DATA  TAPE  FORM  FOR  THERM0  DATA  OUTPUT 
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SECTION  V 


THERAfa  DATA  PROGRAM 


The  program  (Appendix  II)  to  fit  the  tabular  JANAF  thermochemical 
data  to  the  Lewis  polynomial  is  relatively  simple  and  straightfor  ard. 
JANAF  data  for  Cp,  H,  and  S  ir-  read  in  the  tabular  form  for  two  temper¬ 
ature  ranges  298®  to  1,000°K  and  1,000“  to  6,000“K.  The  format  is  as 
shown  in  Table  XII  for  the  gas  TI02.  The  first  line  consists  of  the 
reference  heat  of  formation  to  298“K  in  calories  per  mole,  the  number 
of  points  in  the  first  temperature  range  and  the  number  of  points  in  the 
second  temperature  range,  followed  by  a  species  identifier.  The  tabular 
temperature  references  for  the  first  interval  are  then  listed  followed 
by  the  appropriate  Cp,  H,  and  S  values.  -This  sequence  is  then  repeated 
for  the  second  temperature  range.  Subsequent  to  data  input,  the  program 
establishes  a  matrix  which  it  solves  by  an  inversion  procedure  to  obtain 
the  required  polynomial  constants.  These  constants  are  then  output  in 
the  form  suitable  for  inclusion  on  the  Lewis  program  thermodynamic  data 
tape. 


The  output  is  as  shown  on  Table  XIII  for  the  Ti02  species  which  is 
also  shown  on  the  input.  The  coefficients  required  on  the  data  card  for 
the  first  temperature  range  are  shown  on  the  first  line  followed  by  a 
comparison  between  the  tabular  data  apd  the  generated  polynomials  for 
the  three  functions  of  interest.  It  can  be  seen  that  the  data  fit  is 
quite  good.  These  data  are  followed  by  the  required  constants  for  the 
second  temperature  range  and  a  similar  comparison  of  polynomial  and 
tabular  results . 

These  results  are  read  onto  the  data  tape  in  the  form  shown  on 
Table  XIV.  There  are  four  cards  per  species.  The  first  line  consists 
of  the  empirical  formula  in  columns  1  to  12.  A  gaseous  species  has  no 
subscript,  a  liquid,  has  subscript  L,  and  a  solid,  has  subscript  S. 
Columns  19  to  21  contain  a  code  for  the  data  source,  in  this  case,  a  J 
for  the  JANAF  tables .  Columns  22  to  24  are  the  month  and  year  of  poly¬ 
nomial  computation.  The  first  through  the  fourth  atoms  in  the  species 
are  inserted  in  columns  25  and  26,  30  and  31,  35  and  36,  40  and  41,  while 
the  respective  number  of  atoms  are  placed  in  columns  27  to  29,  32  to  34, 
37  to  39,  and  42  to  44.  Column  45  is  used  for  the  phase  of  the  species 
G,  L,  or  S.  Columns  46  through  55  have  the  lower  temperature  bound  at 
which  the  species  exists  .vith  the  upper  temperature  bound  being  placed 
in  columns  56  through  65. 

Cards  2,  3,  and  4  contain  the  14  polynomial  constants,  five  each  on 
cards  2  and  3  and  four  on  card  4. 
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